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ABSTRACT 

There is mounting evidence in solar wind observations and in numerical simulations that kinetic 
and magnetic energies are not in equipartition in magnetoliydrodynamic turbulence. The origin of 
their mismatch, the residual energy = ~ Eh, is not well understood. In the present work 
this effect is studied analytically in the regime of weak magnetohydrodynamic turbulence. We find 
that residual energy is spontaneously generated by turbulent dynamics, and it has a negative sign, in 
good agreement with the observations. We obtain that the residual energy condenses around fcy = 
with its fc||-spectrum broadening linearly with k±, where fcy and fcj^ are the wavenumbers parallel 
and perpendicular to the background magnetic field, and the field-perpendicular spectrum of the 
residual energy has the scaling Er{k±) oc fcj^ in the inertial interval. These results are found to be 
in agreement with numerical simulations. We propose that residual energy plays a fundamental role 
in Alfvenic turbulence and it should be taken into account for correct interpretation of observational 
and numerical data. 



1. INTRODUCTION 

Turbulence plays an important role in a variety of as- 
trophysical phenomena, including amplification of mag- 
netic fields in planetary and stellar interiors, magneti- 
zation and angular momentum transport in accretion 
disks, scattering of cosmic rays, formation of small- 
scale density structures in the interstellar medium, heat- 
ing of the solar corona and the solar wind. Theoret- 
ical description of magnetized plasma turbulence is a 
complicated task. A valuable guidance for phenomeno- 
logical modeling is pro vided by high-resolution numer- 
ical simulations (e .g., iBrandenburg fc Nordlundl 120111 : 
iKritsuk et al.l 120 lit) . On the analytical side, a signifi- 
cant progress is made in the limit of weak turbulence, 
that is, turbule nce consisting of weakly interacting lin- 
ear w aves (e.g., INewell et all 120011: IGaltier et all I200a 
I2002f ). which provides a testbed for fundamental ideas 
in the theory of turbulence, such as scale invariance, lo- 
cality of interactions, energy cascades, and anisotropy. 
Weak MHD turbulence is also interesting on its own 
as it may play a role in the interstellar medium, in 
the solar corona and the sol ar wind, in planetary and 
stellar magnctosphcrcs (e.g., 'Bhattachariec fc N d 120011: 
Saur et al. 2002: Luo & Melrose 2006: Rappazz o et al.l 
20081: 1013^^^120101: IChandran et al.ll2010l) . 



In the present Letter we address the problem of mis- 
match between kinetic and magnetic energies recently re- 
ported in observations of the solar wind turbulence. Our 
interest is also motivated by a similar mismatch found in 
recent numerical simulations of magnetohydrodynamic 
(MHD) turbulence, suggesting that the phenomenon 
may, in fact, have a fundamental nature rather than re- 
flect a possible non-universality of solar wind turbulence. 
By employing the framework of weak MHD turbulence 
we present the first direct analytical derivation of residual 
energy generation in MHD turbulence. We demonstrate 
that kinetic-magnetic equipartition gets spontaneously 



broken by interacting random Alfven waves, even if it 
is present initially. Our results indicate that this effect 
plays a fundamental role in the energy cascade, and it 
should be taken into account in phenomenological mod- 
eling of MHD turbulence. 

For our analysis we write the incompressible MHD 
equations in terms of the Elsasser variables: 
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where the Elsasser variables are defined as = v ± b. 
V is the fluctuating plasma velocity, b is the fluctuating 
magnetic field normalized by ^/47rpo, = Bo/^/^^rpo is 
the Alfven velocity corresponding to the uniform mag- 
netic field Bq, P = {p/po + ^^/2) includes the plasma 
pressure p and the magnetic pressure, po is the con- 
stant mass density, and we neglected the terms describ- 
ing viscous and resistive dissipation. One observes that 
when zT(x, i) = 0, an arbitrary function z^(x, i) = 
F^{x ± VAt) is an exact solution of ([T]), which repre- 
sents a non-dispersive Alfven wave propagating paral- 
lel or anti-parallel to Bg with the Alfven speed. Non- 
linear interactions are the result of collisions between 
counter-propagating Alfven wave packets. This qual- 
itative picture plays a central role in phenomenologi- 
cal models of MHD turbulence; comprehensive reviews 
of recent analytical and numerica l resul t s can be found 
in. e.g.. (IGaltied [200l ISridhad IMot IMininnil [20Tll: 



INewefl fc Rumpai2011l : iTobias et al 

The ideal MHD system ([T]) conserves the two inde- 
pendent Elsasser energies, £:± = (|z±|2)/4, related to 
the total energy and cross-helicity, E = £'+ + E^ and 
He = E^ — E^ , respectively. In a turbulent state, when 
energy is supplied to the system at large scales, both 
energies E^ cascade toward small scales where they are 
converted into heat by viscosity and resistivity. Theories 
and numerical simulations of MHD turbulence address 
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the Fourier spectra of the energies E^{k) in the inertial 
range of scales, that is, scales much smaller than the forc- 
ing scales and much larger than the dissipation scales. 

In the picture of counter-propagating Alfvcn modes it 
is often assumed that energy spectra are the same for 
both kinetic and magnetic fluctuations, since such en- 
ergy equipartition holds for individual Alfven wave pack- 
ets. However, recent solar wind observations demon- 
strate contradictions with this assumption: the Fourier 
energy spectra of magnetic and velocity fluctuations 
appea r to have different slopes in the inertial inter- 
val ( PodestaeFaplpO?: Sa lem et al]l2009l: iTessein et all 
200a,Chen et al.l2011b;, Wicks et aLlfeoilHLi et al.ll2011l: 



Chen et aLl l2011aD . This is somewhat similar to the 
mismatch in the spectra found in numerical simu l ations 
of strong MHD tur bulence (jMiiller fc GrappinI 120051 : 
iBoldvrev et~alll2011[ ). Most intri guing, however, is the 
fact that significant difference between magnetic and ki- 
netic energies was detected in numerical s tudies of weak 
MHD turbulence (jBoldvrev fc Pered [20091) consisting of 
mostly independent Alfven waves for which equipartition 
between magnetic and kinetic energies holds exactly. 

The mismatch between magnetic and kinetic energies, 
the so-called residual energy Erjk ) = E„{k) — Eh{k) 
was noted in previous works, e.g., (IGrappin et al.l 119831 : 
iZank et all Il996t iMuller fc GrappinI |2005[) . however, it 
has not been fully explored in models of MHD turbu- 
lence, since such models are often based on the conserved 
quantities, E^, E, He, which do not contain information 
about the residual energy. In this work we present the 
analytic derivation of the residual energy in the case of 
weak MHD turbulence. We demonstrate that residual 
energy is spontaneously generated by interacting Alfven 
waves even if it is absent initially, it plays a fundamental 
role in the turbulent dynamics, and it provides a qual- 
itative explanation for the breakdown of magnetic and 
kinetic equipartition observed in the solar wind. 

2. THE ORIGIN OF THE RESIDUAL ENERGY 

We concentrate on turbulence of shear- Alfven waves 
whose polarizations are normal to the guide field, _L 
Bq. We assume that the guide field is strong; its strength 
with respect to the magnetic and velocity fluctuations is 
characterized by a small parameter e ~ z'^/va ^ 1- We 
also assume that the turbulence is weak, that is, the wave 
spectrum in the field-parallel direction is broad enough to 
ensure that the nonlinear terms are small, k^^VA 3> k±z^, 
where k\\ and k± are typical wave numbers character- 
izing the spectral width in the field-parallel and field- 
perpendicular directions and are typical (rms) values 
of the fiuctuations. The Elsasser Fourier energies are 
given by the correlation functions 

(z±(k,t).z±(k',t)) =e2„^q±(fc||,fc_L,<)'5(k + k'), (2) 
while the residual energy is related to 



ergy define d as E^jk) = j^v\er(Je)2Trk± . In previous 
treatments (jGaltier et al.l [2000l l2002f ). the residual en- 
ergy was assumed to be zero as it would be the case for 
independent z+ and z~ waves. In the present study we 
do not make such an assumption, rather, we derive the 
equation for the residual energy. 

The derivation is performed pert urbatively in e, fo l- 
lowing the standard procedure fe.g.. lGaltier et alll2002|) . 
When the nonlinear terms in Eq. ([1]) are neglected, both 
z^(k) oscillate in time and so does the residual energy, 
Q^{k,t) = exp(2ifc||Wyii). When the nonlinear in- 

teraction is taken into account, acquires small correc- 
tions that one finds, to the first nonvanishing order, by 
integrating the equation: 

dtq''{k,t) = -e-2'fc|i"-''*e2 j^^^^j^-h(q)[^-(p) _ ^-(j,)] 

+9-(q)[g+(p) - g+(k)]}5((Z||)<5(k - p - ci)d^pd\.{A) 

where R^^pq = {■nvAl2){k±_ x qx)^(k_L • P±)(ki • 
c\±^ / {k\p\q\) , and we introduce a "smeared" delta- 
function <5(fc||) = (1 — cxp(— 2ifc||t)^t))/(i7rfc||), which ap- 
proaches 5{k\\) + {i / t:)V {1 1 k\\) at large t. Since we are 
interested in spontaneous generation of the residual en- 
ergy, we assume that it is absent initially, therefore, in 
the rhs of Eq. @ we neglected the terms proportional 
to 5°. The derivation of Eq. (U) is standard although 
lengthy, and it will be presented elsewhere. In what fol- 
lows we discuss the main consequences of the obtained 
result. 

3. THE SPECTRUM OF THE RESIDUAL ENERGY 

First, we notice that the integral in the rhs of Eq. Q is 
essentially nonzero. It therefore describes the generation 
of the residual energy by Alfven waves. The important 
conclusion is that the residual energy is generated by 
nonlinear interactions even if it is zero initially. More- 
over, one can directly verify that for all the spectra of 
interest, the integral in the rhs of Eq. (|4]) is positive. As 
we shall see below, this means that the generated residual 
energy is negative, that is, the magnetic energy exceeds 
the kinetic oneQ 

A simple evaluation of the integral in Eq. @ can be 
performed for the case of balanced MHD turbulence, that 
is, for 9+ = q~ . In this case, the El sasser energy spectra 
have the well known scaling (e.g., iNg fc Bhattacharied 



[1991 IGaltier et al.ll2000l) : q^ 



f{k^\)kj_ , where 



/(fc||) is arbitrary, and we assume that it is smooth at 
least initially. For the reasons that will be clear momen- 
tarily, we estimate the integral in the rhs of Eq. @ at 
fell = 0: 



Rk^pq{. . .}6{qi\)d{k p qjd^pd^q 



ak I '^VA, 



(5) 



where the numerical factor is a ~ P{0)E Non- vanishing 
(z+(k,t) • z"(k',t)) ^ e^v\Q°{ki\,k±;t)S{k + k'). (3) of this term at fc|| = implies that Eq. Q has a solution 



In these formulas we assume spatial homogeneity and av- 
erage over a statistical ensemble. The Elsasser energies 
are real, while the residual correlation function satisfies 
Q'^{k\i,k±)* = Q°(— fc||,fcx). Its real part is then the 
residual energy er{k) = ReQ''. In practical applications 
one also uses the phase- volume compensated residual en- 



^ This result is a statistical effect of many independent collisions 
of Alfven waves. It does not imply that a particular collision of two 
Alfven wave packets necessarily produces negative residual energy. 

A more precise numerical integration gives a « 0.0657r/2(0), 
however, our discussion of scaling properties of the residual energy 
does not require the knowledge of numerical coefficients with such 
precision. 



diverging in time, which is unphysical. The reason is that 
when g° becomes large one cannot neglect the nonlinear 
terms containing in Eq. ([4]). As we shall see, in order 
to obtain the convergent result, it is enough to retain the 
terms linear in , which are proportional to q^q^ . Those 
terms describe relaxation of the residual energy due to 
its interaction with the Elsasser energies. Their structure 
can be found from a simple estimate, if one notices that 
these terms should have the same dimensional structure 
as the rhs term in Eq. (j4]), if one replaces one of q^ by q^ . 
The corresponding relaxation term can thus be added to 
the rhs of Eq. (j4|) in the form 75", with the relaxation 
rate 7 ~ —e^VAk^q^^k). The estimate of the relaxation 
rate at k\\ =0 gives: 

^{Q,k^)^-Pe'k^VA, (6) 

where /3 - /(O)0 

Collecting the results ([5]) and ([6]), we rewrite Eq. (j4]) 
in the form: 

dtq°ik,t) = -e^pk^VAq°ik,t) - e^akl^VAe-^'''^^''^* (7) 
which has to be solved together with the initial condition 



0. The solution is given by 



q"{k,t) 



ae^k~^ 



/3e^kA 



2ih 



-2ik^VAt 



and the residual energy is found as 



alSe^kl^ 



d^e^kl 



-ikl 



(8) 



(9) 



The residual energy can be re- written as 

e^(fc||,fcx) = -a7re2fcJ^A(2fc||), (10) 

where the A function has support inside a narrow wedge- 
shaped domain: 

fc|| < /3e^k±/2. (11) 

According to the result \imf^oe/{x'^ + e^) = tt6{x), 
the function A would turn into the (5-function in the 
limit e — > 0. This explains why the region around = 
was relevant in our estimates ([5]) and (|6]). We will refer 
to such a spectrum as "condensate" . It is also convenient 
to define the field-perpendicular spectrum of the residual 
energy. 



Er{k±) 



'^u^er(fc|| , fcj_)27rfc_Ldfc|| = — Q!7r^e^v^fc](i2) 



Expressions (| WT^ are the main results of the present 
work. 

4. NUMERICAL RESULTS 

In this section we test our analytic predictions in 
numerical simulations. In the presence of a strong 
guide field, the universal properties of MHD turbulence 
can be studied in the framework of Reduced MHD, 
which can be effectively used in numerica l simulations of 
both strong and weak turbulence, e.g., (jPmitruk et al.l 

^ It is easy to see that 7 coincides with the inverse time of 
nonUnear spectral energy transfer in weak MH D turbulence, (e.g., 
INg fc BhattachariedllQQg: iGaltier et al.|[2000h . 



2OOI lOughton et all [200I IGaltier fc ChandranI [20061: 
Perez fc BoldvrevI l2008t iRappazzo et al.ll2010D . In our 

numerical set-up, the z"*" and z~ fields are independently 
driven by Gaussian distributed random, short-time cor- 
related forces with zero mean and prescribed variances, 
so that no residual energy is supplied by the driving rou- 
tine. The forces are applied at fc|| = 1, . . . , 16 & fc^ = 1,2 
in the Fourier space. The box is elongated in the field- 
parallel direction according to L\\/L±^ = Bo/brms- We 
employ a fully dealiased pscudo-spcctral method; the 
run has the resolution of 512'^, the Reynolds number 
Re = 2400, and it is averaged over 41 large-scale eddy 
turnover times. We choose viscosity equal to magnetic 
diffusivity. The details on the code and the n umerical 
set-up can be found in (jPerez fc BoldvrevllMoh . 

We find that the residual-energy condensate is indeed 
generated in weak balanced MHD turbulence, as shown 
in Fig. ([1]). The condensate is negative, and its broad- 
ening with k± shown in Fig. ^ is consistent with the 
linear law predicted in Eq. pT|). The scaling behavior 
of the field-perpendicular spectrum of the residual en- 
ergy is also consistent with the analytic prediction ((T^ . 
see Fig. When we repeated the computation for the 
imbalanced case, that is, the case when g"*" 7^ q~ (not 
presented here), we did not detect essential changes in 
the structure and scaling of the condensate. 








Figure 1. The field- parallel spectra of magnetic (solid line) and 
kinetic (dashed line) energies at k± = 6 for weak balanced MHD 
turbulence. The magnetic energy exceeds the kinetic one at small 
Note the logarithmic scale on the vertical axis. 



5. DISCUSSION 

In this Letter wc demonstrated that residual energy 
is spontaneously generated by nonlinearly interacting 
Alfvcn waves, and it accumulates in a narrow region 
of phase space pT|) . This explains the phenomenon 
of 'conden sate' of the residual ene rgy numerically ob- 
served in ([Boldvrev k PerezI [20091) . Although the to- 
tal (phase-space integrated) residual energy is small, 
its presence is crucial for the turbulent dynamics, for 
the following reason. Inside region ([TT]) . that is, for 
fc|| « the residual energy is comparable to the Elsasser 
energies and it follows the same scaling, er{0,k±) ~ 
q'^{0,k±) cx kj^. Since the turbulent cascade of Alfvcn 
waves crucially depends on the modes with (e.g. , 
Montgomerv fc Turned [T98ll: IS ridhar fc Coldreichll 19941: 



Ng fc Bhattacharieelll996l : IGaltier et al.ll200a 120021) . the 

residual energy crucially affects the dynamics of the El- 
sasser modes. 

More formally, this is expressed in the fact that one 
of the Elsasser energies enters the kinetic equations in 
the combination g^(A:|| , fcj_)5(fc||), where (5(fc||) are con- 
centrated in the region pT|). see, e.g., Eq. (|4]). In previous 
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Figure 2. Upper panel: The (normalized) residual energy 
Er{k^^, kj_)/ Er{0, k±) plotted as a function of fcy and k± for 
the weak balanced MHD turbulence. Lower panel: The field- 
perpendicular spectrum of the residual energy for weak balanced 
MHD turbulence compensated by k±. At small k± the deviations 
are strong due to proximity to the forcing region (not shown). 



treat ments of weak MHD turbulence (jGaltier et al.|[2000l 
I2002D . it was assumed that q'^{k\\,kj_) are smooth func- 
tions of fc|| , that is, their dynamically important compo- 
nents with fc|| = have the same fcx-scaling as the com- 
ponents with ^ 0. Such an assumption would be self- 
consistent if the residual energy were absent. However, 
as we have demonstrated, the residual energy is spon- 
taneously generated by interacting Alfven waves. One 
can argue that the presence of the residual energy mod- 
ifies the spectra of the Elsasser energies that, together 
with their smooth parts, = /^(fc||)fcj'^, now acquire 



e''A*(fc|| )fc_L , where 



their own singular parts Sq 

A='=(fc||) are concentrated in region (|lll) . When multi 

plied by 5{k\\) and integrated over both the smooth 
and the singular parts provide comparable contributions 
to the integrals. 

We established that the spontaneously generated resid- 
ual energy is always negative. Although this result is 
obtained in the framework of weak MHD turbulence, 
it provides the first analytic explanation for the obser- 
vational and numerical findings that magnetic energy 
exceeds kinetic e nergy in the inertial interval of MHD 
turbulence, (e.g., iMiiller fc Grappinll2005l : [Podesta et al.l 
[2007t IBoldvrev fc Perezll2'009l) ! 

We derived that the residual energy has the field- 
perpendicular spectrum Er{k±) oc kj^ (|12p . This rela- 
tively shallow spectrum holds in the inertial interval and 
breaks down at sufficiently large k± when the nonlinear 
broadening of the residual energy spectrum in the field- 
parallel direction (fTT|) becomes comparable to the width 
of the field-parallel energy spectra of the Alfven waves, 
g^(fc||). It is easy to see, however, that this is precisely 
the scale beyond which the weak interaction approxima- 
tion breaks down, and the turbulence becomes strong. In 
our future work we will extend our analysis of residual 
energy to the case of strong MHD turbulence. 
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